The scintillation and ionization yield of liquid xenon for nuclear recoils 
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Abstract 

XENON 10 is an experiment designed to directly detect particle dark matter. It is a dual phase (liquid/gas) xenon 
time-projection chamber with 3D position imaging. Particle interactions generate a primary scintillation signal (S 1) 
and ionization signal (S2), which are both functions of the deposited recoil energy and the incident particle type. We 
present a new precision measurement of the relative scintillation yield -£ e // and the absolute ionization yield <3 V , for 
nuclear recoils in xenon. A dark matter particle is expected to deposit energy by scattering from a xenon nucleus. 
Knowledge of •£«// is therefore crucial for establishing the energy threshold of the experiment; this in turn determines 
the sensitivity to particle dark matter. Our -£«// measurement is in agreement with recent theoretical predictions above 
15 keV nuclear recoil energy, and the energy threshold of the measurement is ~ 4 keV. A knowledge of the ionization 
yield Q v is necessary to establish the trigger threshold of the experiment. The ionization yield Q y is measured in two 
ways, both in agreement with previous measurements and with a factor of 10 lower energy threshold. 

Key words: Dark Matter, Liquid xenon, Time-projection chamber, Scintillation quenching, Nuclear recoil 
PACS: 61.25.Bi, 29.40.Mc, 28.20.Cz, 95.35.+d 



1. Introduction 

1.1. Expected signal in XENON10 

There is abundant evidence for a significant cold dark 
matter (CDM) component in the universe (JIHHD, anc ^ 
perhaps the best-motivated candidate is the lightest neu- 
tralino from super-symmetric (SUSY) extensions to the 
Standard Model |4|. A neutralino is expected to be 
non-relativistic and stable, and is more generally clas- 
sified as a weakly interacting massive particle (WIMP). 



"XENON 10 Collaboration 
'Corresponding author 
Email address: pfsOhet . brown . edu (P. Sorensen) 

Preprint submitted to Elsevier 



The open question of the expected mass and cross- 
section of WIMPs is being addressed by numerous di- 
rect and indirect detection experiments ]5]|6]|7), includ- 
ing XENON 10. 

The XENON 10 detector is a liquid xenon time- 
projection chamber. It is designed to directly detect 
galactic WIMPs which scatter elastically from xenon 
nuclei. With velocities of order 10~ 3 c, the recoil en- 
ergy spectra WIMPs with a mass 100 GeV/c 2 incident 
on xenon is predicted to be a featureless exponential 
falling one decade every 30 keV nuclear recoil energy 
(keVr). A particle interaction in liquid xenon creates 
both excited and ionized xenon atoms |8|, which re- 
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act with the surrounding xenon atoms to form excimers. 
The excimers relax on a scale of 10~ 8 s with the release 
of scintillation photons. This prompt scintillation light 
is detected and referred to as the S 1 signal. 

An external electric field (Ed - 0.73 kV/cm) across 
the liquid xenon target causes a portion of the ion- 
ized electrons to be drifted away from an interaction 
site. The electrons are extracted into the gas phase by 
a stronger electric field (~ 10 kV/cm) and accelerated 
through a few mm of xenon gas, creating a secondary 
scintillation signal. This scintillation light is propor- 
tional to the number of ionized electrons and is referred 
to as S2. The amplification during proportional scintil- 
lation makes the recoil energy threshold for S2 lower 
than the threshold for S 1. XENON 10 discriminates be- 
tween electron recoil background and the expected nu- 
clear recoil WIMP signal via the distinct ratio of ioniza- 
tion (52) to scintillation (SI) for each type of interac- 
tion. 



1.2. Importance of these measurements 

The energy threshold of XENON 10 is determined by 
its total light collection efficiency for primary scintil- 
lation photons (SI), and by the effective scintillation 
yield of nuclear recoils (-£ e //)- Because of the expo- 
nential slope of the expected signal, the detector energy 
threshold bears significantly on the ultimate sensitivity 
of XENON10. The sensitivity of XENON10 to spin- 
independent interactions [9| and spin-dependent inter- 
actions 1 10] are reported in separate letters, based on a 
constant -C e // = 0.19. Several groups have measured 
-£?// using tagged neutron scattering, with a range of 

results rrnrrei . 



In Sec. 3.1 we present an alternative method to mea- 
sure £eff. We clearly establish the energy dependence 
of X.eff in the range 4-100 keVr. The uncertainty is 
substantially reduced compared with previous measure- 
ments. In Sec. 13.21 - Sec. |3.5| we make a careful assess- 
ment of possible systematic and statistical uncertainties 
affecting our measurement. The effect on the dark mat- 
ter sensitivity of XENON10 is discussed in Sec. B] 

In Sec. |4] we report a measurement of the absolute 
ionization yield (Q y ) of nuclear recoils in liquid xenon. 
Our results are in agreement with previous measure- 
ments, above 25 keVr Ifl4l . To our knowledge, this is 
the first measurement of (2, below 25 keVr. In Sec. 4.1 



we also present a new method to determine the absolute 
ionization yield. This method provides a cross-check on 
our measurement of X e //. 



2. Experimental Apparatus 

2.1. XENON 10 detector and neutron calibration 

XENON10 is a position-sensitive time-projection 
chamber. Two arrays of UV-sensitive photomultiplier 
tubes (PMT) detect the SI and S2 signals. The 
XENON 10 instrument, including design, energy cali- 
bration and position-dependent corrections, is described 
in detail in 1151 . The performance of the 3D position 
reconstruction is described in [16). The XENON10 de- 
tector is shown schematically in Fig. [T] 
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Figure 1: A side view cut-away schematic of the XENON 10 detector 
as rendered by the GEANT4 simulation. The liquid xenon target is 
15 cm in height with a 20 cm diameter. Ancillary systems, cabling, 
shielding etc. are omitted for clarity. The detector was completely 
enclosed by 20 cm Pb outside of 20 cm polyethylene shielding. 

The XENON10 detector was exposed for 12 hours, 
with a live fraction 0.92, to a 3.7 MBq ±15% AmBe 
source emitting 220 n/s. The neutron rate is based on a 
yield of 6 x 10~ 5 n/Bq |fT7| . The exposure occurred in 
low-background conditions at Laboratori Nazionali del 
Gran Sasso. With 3100 meters water equivalent rock 
overburden, the cosmic muon flux is reduced by about 
10 6 compared with surface conditions [18]. The instru- 
ment was shielded by 20 cm Pb outside of 20 cm high 
density polyethylene. The shield completely enclosed 
the detector. It reduced the cavern y flux by more than 
10 5 , and the cavern neutron flux by about 10 2 5 iPTS! . 

2.2. Neutron source and Monte Carlo simulation 

The AmBe source was attached to a steel rod and in- 
serted through a 7 mm diameter hole in the shield. It 



was positioned next to the detector cryostat, behind an 
additional 5 cm of Pb shielding. The active target of 
XENON 10 had a 10 cm radius and 15 cm height, with a 
xenon mass of 13.7 kg. The analysis presented here uses 
only nuclear recoils which occurred in a 5.4 kg fiducial 
target, with 8 cm radius and 9.3 cm height. This is the 
same fiducial target used for the blind analysis of the 
WIMP-search data, as reported in [9|, surrounded by a 
minimum of 2 cm of self-shielding xenon. 

Initial neutron energies from the AmBe source ranged 
from below 0.1 MeV to 11 MeV, with a mean at 4.3 
MeV 1 19]. The neutron energy spectrum is known with 
an accuracy of ±3% (per 0.1 MeV bin) for a source of 
this strength (20). A detailed Monte Carlo simulation 
of the nuclear recoil spectrum in the xenon target was 
found to be insensitive to variations on this scale. De- 
spite the features in the AmBe source energy spectrum 
lfT9l . the spectrum of neutron energies as they enter the 
5.4 kg fiducial target is a featureless exponential, falling 
1 decade in 3.5 MeV. 

The source also emitted 148 y/s at 4.4 MeV from the 
de-excitation of the 12 C final state. A 0.06 MeV y with 
a branching ratio of 36% was not relevant due to the 
5 cm of internal Pb shielding. The Monte Carlo simula- 
tion predicted a flat rate of single-scatter y events in the 
fiducial target with E < 100 keV electron recoil equiv- 
alent energy (keVee). The predicted rate was reduced 
by x40 by the 5 cm of internal Pb shielding, and by an 
additional x3 due to the 2 cm of self-shielding xenon. 
Prior to S2/S 1 discrimination, the single scatter y rate 
in the energy range E < 100 keVee was measured to 
be < 2 cts/0.25 keVee over the full exposure. In con- 
trast, the single scatter elastic nuclear recoil rate was 
higher by as much as x400, as shown in Fig. [2] Inelas- 
tic nuclear recoils on xenon cause additional prominent 
y lines at 40 and 80 keVee. These events were easily 
avoided by the S2/S 1 discrimination, as shown in Fig. 

El 

2.3. Data acquisition and analysis 

Event data were recorded at a rate of 6.5 Hz during 
the 12 hour exposure. An 52-sensitive trigger was ob- 
tained from the shaped (r = 1 fis) sum of the 30 central 
(of 47) top-array PMTs. The trigger was verified with 
a known voltage input to have full efficiency for an S 2 
corresponding to 4 electrons, and about 75% acceptance 
for 3 electrons. Analysis was restricted to 52 pulses 
corresponding to at least 8 electrons, which is more 
than 3cr below the typical S2 pulse at SI =3 photo- 
electrons. This is shown in Fig. |4] The S 1 signal was 
found by look-back after events were digitized. The S 1 
peak-finding efficiency depends on the size of individual 



photo-electron pulses from the PMTs, and the require- 
ment of a minimum of 0.35 photo-electrons/PMT in 2 
or more PMTs (n > 2 coincidence) in a 0.3 /-is time win- 
dow. The S 1 peak-finding efficiency was predicted to 
be better than 0.985 at 1 keVee. This is discussed fur- 



ther in Sec. 3.3 and is indicated by the dotted line in 



Fig. [5] A high energy veto was set to avoid digitizing 
y scatters with energies above ~ 150 keVee. The veto 
was about xl.5 above the high-energy tail of the elastic 
nuclear recoil spectrum. 

3. Scintillation Yield 

3.1. Measurement of -C e ff 

The nuclear recoil energy E nr (in units of keVr) was 
determined from the measured S 1 response (in units of 
keVee) according to 



Si 



L,y ' -L.eff S n 



(1) 



where -C e // is the scintillation yield of xenon for nu- 
clear recoils, relative to the zero-field scintillation yield 
for electron recoils at 122 keVee. The use of this 
standard reference energy avoids any systematic er- 
ror from the non-linear scintillation yield of xenon for 
electron recoils. After position-dependent corrections, 
XENON 10 measured a volume-averaged light yield 
L v = 3.0 ± 0.1 (syst) + 0.1 (stat) photo-electons/keVee 
for 122 keVee y events 1151 . Scintillation quenching by 
the external electric field Ed - 0.73 kV/cm was pre- 
viously measured to be S e - 0.54 + 0.01 for electron 
recoils and S n - 0.93 ± 0.04 for nuclear recoils IfTD . 

The Monte Carlo recoil energy spectrum was con- 
verted from keVr to keVee, via Eq. [T] and convolved 
with the measured S 1 energy resolution. .£<.// was then 
varied to give the best agreement between the measured 
and Monte Carlo energy spectra, using the binned max- 
imum likelihood method [22|. The primary scintillation 
spectrum of single-scatter nuclear recoils is shown in 
Fig. [2] (black line, with lcr error bars), along with the 
best-fit Monte Carlo spectrum (red line). The fit has 
X 1 = 401 with d.o.f. = 385, and a p-value of 0.28. 
For comparison, the best-fit spectrum assuming a con- 
stant X. e ff - 0.19 is also shown (blue line, x 2 - 549, 
p < 10~ 6 ). To avoid additional systematic uncertainty, 
bins for which the S 1 peak-finding efficiency was pre- 
dicted to be < 0.99 were not included in the maximum 
likelihood fit, and are indicated in grey. The Monte 
Carlo spectrum for these bins is shown dashed. The 
same spectrum without the SI peak-finding efficiency 
is shown dotted. 
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Figure 2: The primary scintillation spectrum (with l<x errors) of sin- 
gle scatter nuclear recoils from 1 1 hrs live exposure to an AmBe neu- 
tron source, and the result from a detailed Monte Carlo of the ex- 
periment (red line, labeled "best fit"), considering the energy depen- 
dence of £. e ff shown in Fig. [31 The S 1 peak-finding efficiency for 
the bin centered on 1.125 keVee (3.4 photo-electrons) is > 0.99. Bins 
< 1 keVee were not used to obtain the fit, and are shown in grey (data) 
or dashed/dotted red (Monte Carlo). The translation from keVr to 
keVee is shown for the 7 X e // spline points from Fig. [3| (red circles 
along axis). Also shown are the Monte Carlo spectra without the ef- 
fect of the S 1 peak-finding efficiency (red dotted), and for a constant 
Leff = 0.19 (solid blue line). 



The form of X. e /f was modeled by a piece- wise cubic 
spline, which provides optimal freedom with continu- 
ous 1 st and 2 nd derivatives. Spline points were fixed at 
2, 5, 10, 15, 25, 50 and 100 keVr and unconstrained in 
-Leff- Spline endpoints (also unconstrained in X. e ff) at 
0.8 and 500 keVr are not shown; dX«// /d{E nr ) 2 - at 
the endpoints. Using fewer spline points was found to 
over-constrain the result. Conversely, using more spline 
points merely added fluctuations (within the statistical 
error). If the keVr values of the spline points were trans- 
lated in either direction, the result shifted by less than 
±0.005. This effect is included in the statistical uncer- 
tainty as shown in Table [T] and Fig. [3] 

The energy dependence of .£>// is shown in Fig. [3] 
(red lines), along with the spline points (red circles). 
Also shown are data from [ 1 1 1 (triangles), lfT2l (squares) 
and the theoretical prediction of [36 1 above 10 keVr 
(dash-dot line). The keVee equivalent value of each 
spline point is indicated in Fig. [2](red circles, along the 
axis). The lcr (stat) error bars on the best-fit -£ e // spline 
points were estimated by a Monte Carlo method that 
is standard for multi-parameter fits: the bin counts of 
the best-fit recoil spectrum shown in Fig. [2] were al- 
lowed to Poisson-fluctuate. The resulting spectrum was 
then treated as input data, and fit with the same method 
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Figure 3: The energy dependence of -Ceff at £rf = 0-73 kV/cm (red 
curves) obtained from the measurement described in Sec. |3.1| and 
the spline points (red circles) with lcr (statistical) errors. The three 
curves indicate the range of systematic uncertainty. The uncertainties 
are summarized in TablefT] Also shown are data from 1111 (triangles), 
[12] (squares) and the theoretical prediction of |36| above 10 keVr 
(dash-dot). -Leff is shown dashed below ~ 3.5 keVr (1 keVee), to 
indicate the energy regime not used to obtain the best fit. A constant 
JLeff = 0.19 is also shown (solid blue line). 



described above. About 10 4 such "experiments" were 
performed, and the resulting .&.// values at each of 
the spline points were found to be normally distributed. 
The overall event rate normalization between the Monte 
Carlo and data was treated as a free parameter in the 
■£?// minimization. Allowing for the DAQ dead time of 
8% and the event acceptance described in Sec. 



3.2 the 



Monte Carlo predicted an absolute single scatter neu- 
tron rate consistent with the uncertainty on the AmBe 
source strength. 

3.2. Event Acceptance 

Event records of 156 [is were digitized with a 
105 MHz ADC. Valid events were required to consist 
of a single 51 pulse followed by a single S2 pulse. 
The pulses are easily distinguishable by their width and 
shape: SI decays with t — 27 ns [23 1, while 52 is 
roughly Gaussian with FWHM ~ 0.8 yt/s. There was 
no distinguishable S 1 pulse for 7% of events inside an 
8 cm radius (the z-coordinate is indeterminate without 
the S 1). A further 4% of events had an S 1 candidate 
pulse with all its signal concentrated in a single PMT 
Such events fail the n > 2 coincidence requirement. 
Both of these observations were expected due to Pois- 
son fluctuations in the small number of photo-electrons 
resulting from low-energy nuclear recoils. In particular, 
the number of events considered as having no S 1 was 
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Figure 4: The discrimination parameter logio(S2/S 1) for nuclear re- 
coils in xenon. The color scale (z axis) indicates logio(counts). The 
parameter is Gaussian except for a small (< 2%) secondary popula- 
tion below the -3<x contour. The analysis thresholds for scintillation 
light (SI) and charge (S 2) are indicated by dashed lines. The centroid 
of the distribution for y-induced electron recoils (red line) is expected 
to lie slightly above the center of the 40 and 80 keVee inelastic popu- 
lations 1 3 1 ] - y events with E > 80 keVee are not shown. 



predicted to better than 15% below 1 keVee. This can 
be seen in Fig. [2] as the difference between the best-fit 
Monte Carlo (dashed), and the same spectrum without 
the 5 1 peak-finding efficiency applied (dotted). 

Basic quality cuts placed loose constraints on accept- 
able 51 and 52 pulse widths, and the distribution of 
scintillation light between top and bottom PMT arrays. 
Additionally, the first and last 50 samples of each event 
were required to be featureless. These cuts had a com- 
bined acceptance > 99% for genuine scatters in the fidu- 
cial target, as previously reported for the WIMP-search 
results |9|. A ±3cr cut on the discrimination parameter 
y = logio(52/5 1) was applied to select the elastic nu- 
clear recoils, as indicated in Fig. |4] This cut excluded a 
small (< 2%) non-Gaussian population of events below 
the -3<x band in Fig. |4] Two other quality cuts were ap- 
plied to the data. The first ("excess Photon Noise") re- 
quired that the proportion of signal concentrated in the 
primary 52 pulse be > 0.85; pulses failing this cut tend 
to have an excess of single-electron 52 events (or other 
spurious photon noise), which can complicate the deter- 
mination of pulse parameters. The second ("single SI") 
required that only 1 possible 5 1-like pulse be found be- 
fore the 52 pulse; a secondary 51 candidate was re- 
quired to have at least x0.25 as many photo-electrons 
as the primary. The rate of SI -only pulses in the ab- 
sence of the AmBe source was > 100 Hz, with an ex- 
ponentially falling energy spectrum. The spectrum of 



0.95 




0.85 



• QC0 (0.993) 
A excessPhotonNoise (0.974) 
■ singleSI (0.971) 
— Total (0.937) 
S1 peak-finding 



0.5 



1 2 

keVee 



5 10 20 

[3S1 phe/keVee] 



50 



100 



Figure 5: Event acceptance for single-scatter nuclear recoils in the 
fiducial target. Average acceptance 0.5 - 100 keVee is indicated in 
the legend. The cuts are described in Sec. |3.2| The S 1 peak-finding 
efficiency is discussed in Sec.|3.3| 



these random coincidence events leads to the mild in- 
crease in acceptance of the "single 51" cut with en- 
ergy, as shown in Fig. [5] Below about 3 keVee the 
acceptance rises again due to the x0.25 photo-electron 
threshold for declaring multiple S 1 candidates. Each of 
these cuts had an average acceptance > 0.97 for elas- 
tic nuclear recoils in the fiducial target. The cuts are 
not completely orthogonal, and the combined average 
acceptance for nuclear recoils - including basic quality 
cuts - was found to be 0.937. The measured -£«// was 
not significantly changed if any or all of these cuts were 
removed, however, the quality of the fit was worse by 
10-20%. 

Single-electron 52 pulses were found in about 50% 
of elastic nuclear recoil events, not necessarily corre- 
lated with the event. About 4% of events also exhibited 
a small secondary 52 < 200 photo-electrons. A single 
scatter nuclear recoil with 52 ~ 200 photo-electrons 
(about 8 electrons) would most likely fall below the de- 
tector 5 1 threshold, as can be seen from Fig. [4] Sec- 
ondary 52 pulses of this size or smaller were not treated 
as multiple scatters. 

3.3. Signal Identification and S 1 Acceptance 

A typical 0.75 keVee elastic nuclear recoil event with 
51= 2.3 photo-electrons and 52 = 899 photo-electrons 
is shown in Fig. [6] Baseline fluctuations of up to 
0.10+0.03 photo-electrons/sample were set to zero. The 
threshold was set individually for each PMT based on 
RMS noise. The waveform is shown summed across 
all 89 PMTs. Coherent noise pick-up was efficiently 
rejected because photo-electron pulses were checked on 



individual PMTs. The pulse at 100 fis is a typical single- 
electron S 2, consisting of 22 photo-electrons spread out 
over ~ 0.8 yus. The S 1 peak-finding efficiency shown 



18 






- 


16 






899 phe 


14 






- 


CD 








p-12 




I 


- 


CO 










o 


i 

1 


- 


£ 8 


X 

<D 




- 


? 6 


Q. 


\ \ 


- 






/ \ l \ 




4 


- 


0.25 us -» 


- 


2 




I 




, , i 





30 



40 



50 



60 



70 



80 



90 



100 



Figure 6: An elastic nuclear recoil event at 0.75 keVee. The S 1 
signal of 2.3 photo-electrons at 30 /is is also shown inset X 10 in y and 
X 100 in x. The S2 signal consists of 899 photo-electrons (about 37 
ionization electrons). The small pulse at 100 /is is a 22 photo-electron 
S 2 from a single ionization electron. Baseline fluctuations of up to 
0.10 ± 0.03 photo-electrons/sample were set to zero in the analysis. 

in Fig. |5](dotted) was obtained by a detailed simulation 
of the photo-electron counting statistics, including PMT 



response and sampling effects, as described in Sec. 2.3 



It predicted a flat 100% efficiency above 2 keVee. The 
low-energy Compton-scatter background observed in y 
calibration data was found to be flat in the energy range 
2-20 keVee. 

Under the assumption that this background should re- 
main flat in the range 0.2 - 2 keVee, the predicted S 1 
peak-finding efficiency and the observed roll-off in the y 
calibration data are consistent at the lcr level. In order to 
estimate the possible systematic effect of the S 1 peak- 
finding efficiency on the X. e ff result, an S 1 efficiency 
curve shifted higher (lower) by 20% was also consid- 
ered. The effect on £ e ff at 2 keVr was a systematic 
decrease (increase) of 0.012. At all other spline points, 
the shift was at most half this amount. This is summa- 
rized in Table [T] As an additional check, we found that 
■Ceff shifted by < 0.005 if the analysis threshold of 
1 keVee was lowered to 0.5 keVee. 

3.4. Input to the Simulation 

The Monte Carlo simulation was performed us- 
ing GEANT4.9.0. By default this software uses the 
ENDF/B-VI nuclear cross section data, in which xenon 
cross sections are based on [24 J. Instead, we used the 
updated ENDF/B-VII xenon nuclear cross section data 



Table 1: -£<>// and the statistical uncertainty are tabulated for each 
spline point shown in Fig. 15] Also shown are the systematic errors 
arising from uncertainty in the S 1 peak-finding efficiency ("S 1 eft'.") 
and the Xe(n,n)Xe elastic scattering cross-section ("<x<,/ data"). These 
systematic uncertainties are discussed in Sec.|3.3|and Sec.|3.4| 



keVr 


£eff 




Uncertainty 








statistica 


Sleff, 


o~ e i data 


2 


0.160 


±0.014 


±0.012 


±0.009 


5 


0.156 


±0.011 


±0.006 


±0.001 


10 


0.162 


±0.012 


±0.005 


±0.002 


15 


0.194 


±0.011 


±0.006 


±0.005 


25 


0.220 


±0.012 


±0.005 


±0.001 


50 


0.237 


±0.009 


±0.004 


±0.001 


100 


0.274 


±0.010 


±0.004 


±0.005 



to obtain the nuclear recoil spectrum, which is based on 
more recent calculations using the Optical Model Po- 
tential (OMP) parameters in 11251 l26l . The uncertainty 
in our Monte Carlo recoil energy spectrum arising from 
uncertainty in the o~ e \ data for Xe(n,n)Xe scattering was 
therefore determined by the OMP theory. 

A comprehensive study of calculated versus mea- 
sured differential neutron cross sections of nuclei from 
A - 19 to A - 209 [29 1 allowed the refinement of the 
OMP theory for nuclei in that mass range. A related 
study [ 28 ] concluded that for nuclei where experimental 
elastic scattering data is lacking, a ±3% uncertainty in 
the OMP potential well depth parameter could be used 
to determine a conservative estimate of the uncertainty 
in the elastic cross sections. We used the EMPIRE soft- 
ware package [ 30 1 along with the best OMP model for 
xenon to calculate d'cr^i/dEdQ for all relevant neutron 
energies. The cases -3% and +3% in the OMP potential 
well depth parameter were considered, and new nuclear 
recoil energy spectra were obtained. An average change 
of < 2.5% in the nuclear recoil spectra was found, in the 
energy range 1-100 keVr. 

The X. e ff analysis was repeated for the two new 
spectra. The resulting measured energy dependence of 
JLeff is shown in Fig. [5] (solid red lines) along with 
the result using the default Monte Carlo energy spec- 
trum. This systematic uncertainty in £ e ff due to the 
input elastic scattering cross section data was found to 
be smaller than the statistical uncertainty. Additionally, 
the older ENDF/B-VI cross section data resulted in an 
energy dependence of -C e // that is consistent with the 
result shown in Fig. |3]for 2 < E„ r < 50 keVr, and about 
2a lower at lOOkeVr. 

We also note that the only available experimental 
measurements of the total cross section for neutrons on 



xenon are from [27|, which are consistent with the OMP 
theory used in this paper. However, this experimental 
data provides weak overall constraints on the uncertain- 
ties in the double differential elastic cross section. 

3.5. S 1 resolution 

The S 1 resolution for nuclear recoils is dominated 
by binomial fluctuations in the collection of scintilla- 
tion photons. It is further impacted by the intrinsic 
single photo-electron resolution of the PMTs, which 
had an average cr/fj. = 0.58 ± 0.05. From this, the 
SI resolution was found to be xl.16 larger than that 
expected solely from the binomial fluctuations, in the 
range E nr < 300 keVr PT1 . This result was deter- 
mined from the width of the log w (S2/S 1) band. The 
total light detection efficiency for S 1 was measured to 
be 8 - 11%. Considering the uncertainty in the PMT 
single photo-electron response, and taking a conserva- 
tive estimate of ±20% uncertainty in the S 1 light detec- 
tion, the uncertainty in the S 1 resolution was calculated 
to be < ±2% [ 3 1 1 . If the resolution were better (worse) 
by this amount, the X. e ff curve in Fig. [3] would translate 
higher (lower) by less than 1/4 the statistical error. This 
effect was included in the statistical uncertainty. 

3.6. Behavior of -C e ff 

The expectation that -C e // should continue to fall with 
decreasing energy is due to the Lindhard model |32| for 
the partitioning of recoil energy into atomic motion and 
electronic excitation (so-called nuclear quenching). In 
Ge the agreement between experiment and the Lindhard 
prediction is very good, above 10 keVr 031 . Below 
10 keVr measured values are as much as x2 above the 
Lindhard prediction |34|. 

In Fig. [7] we show the shape of the simulated single 
scatter nuclear recoil spectrum if X e // were to continue 
to drop monotonically below 10 keVr. This was done 
by taking the best-fit -G// shown in Fig. [3] and forcing 
the spline point at 2 keVr to have X. e ff = 0.10. The 
prediction of about 20% fewer events per bin for all bins 
1-3 keVee is more than 5<x inconsistent with the data 
(X 2 > 700), showing that a decreasing £ e // below 10 
keVr is strongly disfavored by our data. 

The Lindhard model applies to the total detectable en- 
ergy transferred to electronic excitation. The Lindhard 
prediction of X e // in xenon is systematically about 0.03 
higher than our result in the energy range 10 - 100 keVr 
(see for example [35]). This is not surprising since 
scintillation light is generated by electrons which re- 
combine. A theoretical model of biexcitonic electronic 
quenching has been proposed [36]. This prediction for 
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Figure 7: A scenario in which JUff continues to fall below 10 keVr, 
to 0.10 at 2 keVr (as shown inset, blue line) is inconsistent with our 
data. The best-fit Monte Carlo (red line) and energy dependence of 
JLeff (inset, red line) are as shown in Fig.pland Fig. [5] 



the total quenching of xenon scintillation is lower than 
the Lindhard prediction (when applied directly to scin- 
tillation), and is shown in Fig. [3](dash-dot line). 

It has also been noted [ 35 37 1 that for nuclear recoils 
below ~ 10 keVr in xenon, the underlying Thomas- 
Fermi interaction potential may be a poor approxima- 
tion, and the Lindhard model for nuclear quenching may 
not be valid. This is roughly where our -£ c // result 
ceases to fall with decreasing energy, and may indicate 
the point at which the electronic stopping becomes com- 
parable to the nuclear stopping 071 . 

4. Ionization Yield 

The ionization yield Q y is uniquely determined from 
the 51 and S2 data, once X. e ff is known. This is 
shown in Fig. [8] (blue stars) for each of the -£ e // spline 
points from Table [T] This result, for E,j = 0.73 kV/cm, 
agrees with previous measurements [ 14] above 20 keVr. 
Data from lfT4l are shown in Fig. [8] as filled and 
open circles/squares, with the uncertainty omitted for 
clarity. The filled and open data points indicate mea- 
surements from 2 different detectors, while the cir- 
cles (squares) represent measurements at Ed - 0.10 
{Ed - 2.0) kV/cm. The modest variation of Q y with 
applied electric field is expected since -£«// (and by in- 
ference Qy) has been measured [ 1 1 J to have only a weak 
(< 10%) dependence on the applied electric field. 

4.1. Multiple Scatter Method 

A single ionization electron extracted to the gas was 
found to correspond to an S 2 signal of 24 ± 7 photo- 
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Figure 8: The ionization yield Q y (blue stars) of xenon for nuclear re- 
coils at Ej = 0.73 kV/cm, as determined from the energy dependence 
of JLeff (shown in Fig. [51. <2 V as determined from the multiple scatter 
method described in Sec. Hlis also shown (red x), along with data 
from 1 14 1 (filled and open circles/squares, uncertainty not shown). 
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Figure 9: (left) The number of single, double, etc. scatters in the fidu- 
cial target with S 2 (per vertex) > 1600 photo-electrons (black circles) 
is compared with the Monte Carlo prediction (red line). Only scatters 
with an energy deposition per vertex E > E,i„- were counted. The 
comparison is a sensitive function of the Monte Carlo energy thresh- 
old. The best fit corresponds to E, hr = 14.2 ± 0.7 keVr. (right) The^ 2 
distribution (red circles) for the comparison shown at left as a func- 
tion of the recoil energy threshold imposed in the Monte Carlo. A 
2"''-order polynomial fit (black line) was used to extract the minimum 
and the l<x uncertainty (indicated by dashed lines). 



electrons ED . 52 pulses were clearly identifiable over 
> 5 orders of magnitude in size, and as many as 4 scat- 
ters in the fiducial target were observed for some elastic 
nuclear recoil events. 

The multiplicity (m = 1...3) of an elastic nuclear re- 
coil event is defined here as the number of scatters with 
recoil energy above a given threshold. The analysis in 
Sec. [3] was restricted to events with m = 1 and an 52 
threshold of 200 photo-electrons, or 8 ionization elec- 
trons. The measured multiplicity is a sensitive function 
of the specified recoil energy threshold per scatter. A 
specific 5 2 threshold (5 2,/, r ) in the data was correlated 
with a keVr equivalent threshold (E t hr) from the Monte 
Carlo, by comparing the predicted and experimentally 
measured event multiplicities. Due to the typical 0.8 /is 
FWHM of 52 pulses, scatters with dz < 3 mm were not 
easily resolvable, regardless of their separation in (x, y). 
This criteria was duplicated in the Monte Carlo analysis. 

An example of the multiplicity comparison is shown 
in Fig. [9](left)for52,/„. = 1600 photo-electrons. In each 
case, only events with all scatters in the 5.4 kg fiducial 
target were counted. The best fit between the number of 
scatters for m - 1...3 in the data and the Monte Carlo 
was obtained for £,/„ = 14.2 ± 0.7 keVr. The overall 
normalization and the recoil energy threshold were free 
parameters, so the minimum x 2 — 2 with d.o.f. - 1 
gives a p- value of 0. 16. The x 2 distribution for the data 
versus Monte Carlo comparison of scatter multiplicity 
is shown as a function of E,i ir in Fig. |9](right). 

The example shown in Fig. ^corresponds to the data 
point (red x) at 14.2 keVr in Fig. [8] The ionization 



yield was calculated as Q y = 1600/24/14.2 = 4.7 ion- 
ization electrons/keVr. This value was then corrected 
to account for the measured S2 resolution, since E t hr 
from the Monte Carlo has perfect resolution. The cor- 
rection increases GL by an average of 11%. The error 
bars include the uncertainty in the mean of the single- 
electron S2 peak, and the uncertainty from each fit (as 
shown in Fig. m right). The^ 2 distributions for all data 
points shown in Fig. [8] exhibited a clear minimum, a 
X 2 /d.o.f. ~ 1 and were well-fit by a 2" i, -order polyno- 
mial. 



4.2. Discussion 

The major benefit of the multiple scatter method de- 
scribed in Sec. Hlis that it provides an independent check 
on the £ e ff result shown in Fig. [3] Above 5 keVr, 
the multiple scatter method (Fig. [8] red x) predicts 
a Q y that is consistent with the values inferred from 
the best-fit -L e ff- Below ~ 5 keVr, the multiple scat- 
ter method predicts a constant Q Y ~ 6 ionization elec- 
trons/keVr (not shown in Fig. IS). In this regime, the 
result of the method is affected by the light collection 
for primary scintillation (SI). This somewhat unex- 
pected result can be understood from Fig. |4] Lines 
of constant 52 lie parallel to the line corresponding 
to 52 = 200 photo-electrons (8 ionization electrons), 
which is shown dashed in the lower left. For example, 
at E t hr - 2 keVr, 52,;, r « 300 photo-electrons and a sub- 
stantial number of events used in the comparison have 
E < 1 keVee. This implies that an additional number of 
events were lost (due to the 5 1 signal acceptance). The 



multiplicity of these events is not known, and this likely 
affects the determination of Q v . 



5. Summary 

The energy threshold of a liquid xenon time- 
projection chamber using proportional scintillation for 
S 2 is presently determined by its sensitivity to primary 
scintillation (SI) and by -£«//. We have presented the 
energy dependence of X. e ff based on an analysis of nu- 
clear recoils in the XENON10 experiment. XENON10 
recently reported a WIMP-nucleon (elastic scattering) 
exclusion limit of a = 4.5 x 1(T 44 (8.8 x 1(T 44 ) cm 2 
for a WIMP mass of 30 (100) GeV/c 2 assuming a 
constant Xe/y = 0.19. A conservative upper bound of 
o- = 5.2 x 10~ 44 (10.4 x 10~ 44 ) cm 2 was quoted in |g] 
based on a preliminary version of this work. Consider- 
ing the energy dependence of JUff shown in Fig. [3] the 
XENON 10 WIMP-nucleon exclusion limit was found 
to be o- = 5.1 x 10~ 44 (10.2 x 10~ 44 ) cm 2 for a WIMP 
massof30(100)GeV/c 2 . 

We also presented the absolute ionization yield of 
xenon for nuclear recoils (<3 V ), as a function of recoil 
energy. For XENON 10, this clearly shows that the trig- 
ger threshold for nuclear recoils is as low as ~ 1 keVr. 
The analysis threshold is presently limited by the light 
collection for primary scintillation (S 1). The technique 
used to determine -£<.// for primary scintillation (SI) 
should also be applicable to the secondary scintillation 
(S2); a future work will explore this possibility. 
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